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ABSTRACT 

We present calculations of non-LTE, relativistic accretion disk models applicable to the high/soft 
state of black hole X-ray binaries. We include the effects of thermal Comptonization and bound- free 
and free-free opacities of all abundant ion species. Taking into account the relativistic propagation 
of photons from the local disk surface to an observer at infinity, we present spectra calculated for 
a variety of accretion rates, black hole spin parameters, disk inclinations, and stress prescriptions. 
We also consider nonzero inner torques on the disk, and explore different vertical dissipation profiles, 
including some which are motivated by recent radiation MHD simulations of magnctorotational turbu- 
lence. Bound-free metal opacity generally produces significantly less spectral hardening than previous 
models which only considered Compton scattering and free-free opacity. It also tends to keep the 
effective photosphere near the surface, resulting in spectra which are remarkably independent of the 
stress prescription and vertical dissipation profile, provided little dissipation occurs above the effective 
photosphere. This robustness is due to the fact that the disk scale height is largely independent of the 
stress and dissipation when the disk is radiation pressure supported and electron scattering dominated 
throughout its interior. Our model spectra may be affected by continuing theoretical uncertainties 
that we do not account for, and we discuss these uncertainties. We provide detailed comparisons 
between our models and the widely used multicolor disk model. Frequency dependent discrepancies 
exist that may affect the parameters of other spectral components when this simpler disk model is 
used to fit modern X-ray data. For a given source, our models predict that the luminosity in the 
high/soft state should approximately scale with the fourth power of the empirically inferred maxi- 
mum temperature, but with a slight hardening at high luminosities. This is in good agreement with 
observations. 

Subject headings: accretion, accretion disks — black hole physics — radiative transfer — X- 
rays: binaries 



1. INTRODUCTION 

The high state spectral energy distribution (SED) of 
black hole X-ray binaries is dominated by a soft ther- 
mal component generally believed to be emission from 
an optically thick accretion flow. It is often hypothe- 
sized that the geometry of this flow is a radiatively effi- 
cient, geometrically thin disk (Shakura & Sunyaev 1973, 
hereafter SS73) or at higher luminosities a less radia- 
tively efficient slim disk (Abramowicz et al. 1998). Most 
sophisticated spectral modeling of black hole accretion 
disks has focused on the case of supermassive black holes. 
However, models that compute vertical structure and ra- 
diative transfer have been constructed for Galactic black 
hole candidates assuming thin disk (Shimura & Takahara 
1995, hereafter ST95; Merloni, Fabian & Ross 2000) and 
slim disk (Wang et al. 1999) relations. 

ST95 argue that a simple diluted blackbody provides 
adequate fits to the local specific flux F v for regions of 
the disk which produce most of the X-ray photons. They 
conclude this for disks with sufficiently high luminosity 
over a range of luminosities, black hole mass, and black 
hole spin. The local specific flux is then given by 



Fu — -j4-B u (f co \T e ff) 

/col 



(1) 
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where B v is the Planck function, T c ff is the effective tem- 
perature of the emitting surface, and f co \ is the color cor- 
rection. They find that in many cases most of the disk 
spectrum is well approximated by a multi-temperature 
diluted blackbody with a single value of / co i. Further- 
more, the value of / co j is not a strong function of the disk 
parameters. They attribute these results to saturated 
Comptonization playing a dominant role and producing 
Wien-like spectral profiles. The relative constancy of / co i 
in these model annuli has been used to justify fitting 
high state SED's with the multicolor blackbody model 
{diskbb in XSPEC) of Mitsuda et al. (1984). Although 
this model assumes a simple profile of effective temper- 
ature with radius, T g oc i?~ 3 / 4 , and does not account 
for relativistic transfer effects, it has been used to infer 
inner radii and black hole spins for disks believed to ex- 
tend deep within the gravitational potential of the black 
hole. 

Merloni et al. (2000) caution against assuming a nearly 
constant value for / co i. They construct several constant 
density disk models with a range of accretion rates. They 
assume some fraction of the released gravitational energy 
dissipates directly into an optically thin corona above the 
disk (Svensson & Zdziarski 1994), but they do not ac- 
count for the irradiation of the disk surface by the coronal 
emission. They fit the resulting SED's with the diskbb 
model, apply relativistic corrections to the best fit pa- 
rameters, and use the results to infer an inner radius for 
the disk. They find that the value of / co i needed to ob- 
tain the correct radius varies significantly from model to 
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model. 

Gierliriski & Done (2004; hereafter GD04) investigate 
temperature and luminosity evolution in a sample of 
black hole binaries observed in the high state with the 
Rossi X-ray Timing Explorer (RXTE). They fit their X- 
ray spectra with a model that includes a diskbb com- 
ponent and accounts for Comptonized emission from a 
corona, focusing on observations in which the Comp- 
tonized emission is < 15% of the bolometric luminosity. 
They find that many of the well observed sources are 
consistent with weak evolution of / co i in which / co ] in- 
creases with higher luminosity or temperature. Further- 
more, they tentatively conclude that such evolution is 
inconsistent with a modified stress prescription in which 
the accretion stress scales with only the gas pressure and 
not the total pressure. 

To our knowledge, all the detailed theoretical models of 
thermal disk spectra of black hole X-ray binaries have so 
far neglected the effects of bound-free opacity of metal 
ions. (This is in marked contrast to models of reflec- 
tion spectra in these systems.) Bound free opacity can 
often dominate free-free opacity at observed X-ray ener- 
gies, reducing the relative importance of electron scatter- 
ing and therefore altering the expected color correction 
factors. This may therefore affect the fits to high/soft 
state spectra and their interpretation. Bound-free opac- 
ity also imparts absorption/emission edges in the spec- 
trum which are potentially observable, particularly with 
the high throughput spectra being produced by modern 
X-ray observatories. 

In addition, significant theoretical advances have been 
made recently in understanding the nature and vertical 
distribution of turbulent dissipation (Turner 2004), as 
well as highlighting the possibility of significant external 
torques on the disk at the innermost stable circular orbit 
(e.g. Gammie 1999; Krolik 1999; Hawley & Krolik 2002). 
These results imply substantial modifications of the basic 
assumptions underlying the disk models upon which all 
spectral calculations have been based thus far. 

In this paper, we construct fully relativistic accre- 
tion disk model atmospheres and examine the integrated 
SED's from these disks. We fully account for bound- 
free opacity by incorporating non-LTE rate equations for 
the ground-state level populations of all abundant metal 
ions. We focus on understanding the sensitivity of these 
SED's to the underlying assumptions in our models. We 
compare our results with previous theoretical investiga- 
tions and recent observations. In sectionElwe review the 
method used to construct our models. We describe our 
standard models in detail in section |3| We report our 
results from altering the stress prescription in section 0] 
adding an inner torque in section [SI and modifying the 
dissipation profile in section [H] In section [7| we compare 
our model SED's with the work of GD04. We provide 
further discussion of all these results and a summary of 
our conclusions in sections and |5| Readers primar- 
ily concerned with the observational implications of our 
models may wish to skip directly to sections |H| and El 

2. METHOD 

We construct models using the methods described in a 
series of papers (Hubeny & Hubeny 1997, 1998, hereafter 
HH98; Hubeny et al. 2000, Hubeny et al. 2001, hereafter 
HBKA). We construct thin disk models by solving the 



fully relativistic one zone disk structure equations in the 
Kerr metric (Novikov & Thome 1973, hereafter NT73; 
Page & Thorne 1974; Riffert & Herold 1985). Next, we 
compute the vertical structure and local spectra for an- 
nuli evenly spaced in log r where r = R/R s , R is the 
Boyer-Linquist radial coordinate, and R g — GM/c 2 is 
the gravitational radius. We compute ten annuli per 
decade in radius and assume the disk extends from the 
radius of marginal stability r ln = r ms to r out ~ 1000. 
Emission from radii larger than r out is assumed to be 
blackbody and contributes little to the integrated flux at 
photon energies > 0.1 keV. 

At each annulus, we use the TLUSTY stellar atmo- 
spheres code (Hubeny & Lanz 1995) to simultaneously 
solve the equations for the vertical structure and angle 
dependent radiative transfer. At X-ray temperatures, 
metal opacities are important, and we incorporate fully 
non-LTE ground state level populations for all ions of H, 
He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni assuming 
solar abundances. Bound- free opacity due to each ion 
is included but bound-bound transitions are neglected. 
We account for Comptonization with an angle-averaged, 
Kompancets treatment of the electron scattering source 
function. Finally, we calculate the integrated disk spec- 
trum seen by an observer at infinity by using a fully 
general-relativistic transfer function (Agol 1997). 

We make several approximations which may be vio- 
lated in nature. The disk models are time steady and 
azimuthally symmetric. All heat is transported verti- 
cally by radiative energy flux. There is no mass loss 
from the disk, so the accretion rate is independent of ra- 
dius. All radiative flux is assumed to originate within the 
disk and its atmosphere; irradiation of the disk surface 
due to an external X-ray source is neglected. Irradiation 
of the disk surface by returning radiation (Cunningham 
1976) is also ignored. The atmosphere calculations are all 
one-dimensional despite evidence that radiation pressure 
supported accretion disks have significant inhomogeneity 
(Turner et al. 2003, Turner 2004). 

For our base model (hereafter the standard model) we 
make additional assumptions which we later relax. First, 
there is no torque on the inner edge of the disk. Second, 
the energy dissipation rate per unit volume e = V • F is 
locally proportional to the density p. Third, we relate the 
vertically averaged stress, r r 0, to the vertically averaged 
total pressure, P using the prescription (SS73) 

f r4> = aP. (2) 

This standard model is parameterized by only four 
quantities: black hole mass M, spin parameter a, the 
stress parameter a, and accretion rate (SS73; NT73). 
The accretion rate is often expressed relative to the Ed- 
dington accretion rate as fa = Mc 2 / XEdd where LEdd = 
1.5 x 10 38 (M/M Q ) ergs -1 is the Eddington luminosity 
for completely ionized H. We prefer to use the quantity 
I = L / L-Edd = rj fa as it relates more directly to the ob- 
served luminosity. Here L is the bolometric luminosity 
of the disk and r)(M, a) is the fraction of the energy at 
infinity which is radiated before crossing r ms . The mass 
is fixed at 10M Q for all models considered here. 

The atmospheres of the individual annuli are entirely 
determined by specifying the midplanc column mass, the 
radiative flux, the local gravity, the composition and the 
vertical dissipation profile e(z). The midplane column 
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mass too is equal to half the surface density of the disk. 
The flux is parameterized by the effective temperature, 
T c ff, of the corresponding blackbody emitter. The verti- 
cal component of the gravitational acceleration g = Q z 
is proportional to the distance z above the midplane and 
the constant of proportionality, Q, is a function only of 
R. These quantities (to , T c s, Q) are uniquely deter- 
mined at a radius r in a one-zone model by a given set 
of disk parameters (M, I, a, a). As noted above, we con- 
sider e(z) oc p(z) for the standard models, but explore 
other prescriptions in section |S| 

Except for the method of calculating m , the stan- 
dard models are equivalent to the supermassive black 
hole models described in HBKA and the reader is re- 
ferred there for further details. HBKA solved a one-zone, 
algebraic relation for too in terms of M, a, a, I, and r. 
This too is then used throughout the stellar atmospheres 
calculation. Due to the simplifying assumptions of the 
one-zone model, there is no guarantee that equation [21 
will hold for the converged solution. The vertically inte- 
grated pressure P can be calculated from the converged 
non-LTE model so the consistency of this prescription 
can be checked after the fact using equation [21 The one- 
zone calculation provides a suprisingly accurate approx- 
imation, but order unity discrepancies do typically arise 
between the final a calculated from the converged non- 
LTE model and the input a.. 

We alter this prescription by using the initial too to 
construct the LTE-grey model (Hubeny 1990) from which 
we calculate an a' with equation [3 If a' is within 10% 
of a we continue the full atmosphere calculation using 
the initial too- Otherwise, we iterate by choosing a 
new too, recomputing the LTE-grey model and a' un- 
til this condition is satisfied. We use Newton-Raphson 
to calculate the change in too for the next iteration: 
SrriQ = dm/da(a — a') with the derivative evaluated ana- 
lytically assuming the one zone scaling holds. The result- 
ing LTE-grey model is then used as the starting point for 
constructing the full non-LTE atmosphere. For most an- 
nuli, this method rapidly converges and provides a final 
a consistent (within ~ 10%) with the input value. How- 
ever, in some of the hottest annuli of the torqued disk 
model, the method fails to converge due to inaccuracies 
in the LTE-grey model. In this case, we default to the 
one-zone calculation for all annuli in the disk, but in all 
other cases we utilize the iterative procedure. 

3. THE STANDARD MODEL 
3.1. The Local Spectrum 

The disks in X-ray binaries with a black hole accretor 
can be much hotter than a disk in a cataclysmic variable 
or quasar. The gas temperatures at the surface of these 
disks are sufficiently high (~ 10 6 K - 10 7 K) that electron 
scattering dominates the opacity at typical photon ener- 
gies and the emitted spectrum may deviate significantly 
from that of a blackbody (SS73). 

We define the effective optical depth as 

rf ee ^ [34 h « s + 4 h )] 1/2 dm' (3) 
Jo 

where K* h — + n]f is sum of the total bound-free and 
free-free absorption opacity, k° s is the electron scatter- 
ing opacity, and to is column mass measured from the 



surface. A closely related quantity is the depth of forma- 
tion 

/>m* 

<= (K™ + K?)dm' (4) 
Jo 

where to* is the frequency dependent column mass at 
which r° ff = 1. All of the disk models considered in this 
paper are effectively thick (i.e. to* < too) at each annu- 
lus for frequencies of interest. In an electron scattering 
dominated atmosphere K* h <C k^ s for typical photon en- 
ergies implying that the photon destruction probability, 
e v = K* h /(>4 h + kI s ) < 1 and r* > 1. Due to the tem- 
perature gradients in the atmosphere, the temperature at 
the depth of formation T* is generally greater than T e g. 
In the absence of Comptonization, the frequency depen- 
dence of e„ and T* alter the spectrum and produce a 
modified blackbody. 

Due to the potentially large number of scatterings 
(n cs w t* 2 3> 1), Comptonization may also be impor- 
tant. To account for its effect on our model spectra we 
always solve the radiative transfer equation including the 
Compton scattering source term. However, we still wish 
to have an indicator of the impact of electron scattering 
on the SED. For a homogeneous electron distribution, 
the y-parameter provides a simple means of characteriz- 
ing the importance of Compton scattering. To gain sim- 
ilar insight, we generalize this prescription by defining a 
frequency dependent y-parameter (HBKA) 



y ;E^fmax[r; 2 ,r;] (5) 

where /cb is Boltzman's constant, T* is the temperature 
at the depth of formation, m e is the electron rest mass, 
and c is the speed of light. 

In Figure 2] we plot the locally emitted SED from a 
single annulus with the frequency dependent quantities 
defined above. The model atmosphere was calculated for 
an annulus at r — 12.6 in a disk with I = 0.1 and a = 0.1 
accreting onto a Schwarzschild black hole. This is near 
where the flux (and thus the effective temperature) of the 
disk peaks. In the top panel (a) we compare the intensity 
at an inclination i of 55° to the vertical with that emitted 
by a diluted blackbody at the same T e g with / co i = 1.56. 
The SED is a modified blackbody with a spectral break 
at the He-like Fe photoionization edge near 8.8 keV. Be- 
low this edge and above ~ 2 keV, the spectrum is well 
approximated by the diluted blackbody. In (b) we plot 
the relative contributions of electron scattering, free-free 
absorption, and bound-free absorption to the total opac- 
ity. The opacities are evaluated at to*. As expected, 
electron scattering dominates except at low energies and 
just blueward of the He-like Fe edge. Bound- free exceeds 
free-free absorption above ~1 keV and therefore should 
not be neglected. Note that a decrease in our frequency 
grid resolution below ~ 6 keV contributes to the breadth 
of the photoionization edges in that range. Panel (c) 
shows how t* depends strongly on the presence of met- 
als. As photon energy increases, the bound-free opacity 
overtakes the free-free opacity, preventing the depth of 
formation from exceeding 10. If we had neglected metal 
opacities, the depth of formation would continue to grow 
as the ratio of free-free to scattering opacity declined. 
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In that case, the photons reaching an observer are cre- 
ated deeper in the atmosphere where the temperatures 
and densities are higher. A direct consequence is seen in 
(d) where y* < 1 for all frequencies of interest. Comp- 
ton scattering may contribute to the broadening of pho- 
toionization edges, but it is not significantly altering the 
continuum. If we had neglected metals, the depth of for- 
mation would have been larger, the number of scatterings 
would have been higher, y* would have been greater than 
unity for hv > 2 keV, and Comptonization would have 
had a larger effect. 

Despite the limited role of Comptonization, the specific 
intensity is still reasonably well approximated by the di- 
luted blackbody spectrum. A modified blackbody caused 
by dominant electron scattering opacity, strong tempera- 
ture gradients, or both is generally not well characterized 
by a single temperature so this result is somewhat sur- 
prising. It seems to be due to the relative constancy of 
r* over the range of photon energies near the peak. This 
means that most of the photons are being created by gas 
at a narrow range of temperatures and that the photon 
destruction probability is only weakly dependent on fre- 
quency. Both of these effects lead to a more blackbody 
like spectral profile but T* is still larger than T c ff. In 
this case, T* « 1.5 — 1.6T c ff near the peak consistent 
with the /coi = 1.56 used for the plot. If we had ne- 
glected bound-free metal opacity and Comptonization, 
the absorption opacity would be more frequency depen- 
dent and the diluted blackbody would be a poorer match 
to the emission. 

Electron scattering and temperature gradients also af- 
fect the angular dependence of the radiation field. As 
shown in Figure[3J the spectra are significantly limb dark- 
ened at frequencies above 0.1 keV and the degree of limb 
darkening is frequency dependent. The curves in Figure 
El show the normalized specific intensity for emission at 
r = 12.6 from an I — 0.1, a = 0.1, and a = disk model. 
The curves are normalized by the specific flux so that 
isotropic emission would have an ordinate of unity. (The 
reason we chose i = 55° in Figure 1 is that the limb dark- 
ened emission matches isotropic emission near this view- 
ing angle.) The curves are plotted for photon energies of 
0.21, 2.2., and 8.4 keV. In this case, the emission near the 
peak (2.2 keV) is well approximated by the results for a 
semi-infinite Thompson scattering atmosphere assuming 
a Rayleigh phase function (Chandrasekhar 1960). It is 
less limb darkened at lower frequencies and more limb 
darkened at higher frequencies due to the temperature 
gradient in the atmosphere. 

3.2. Disk Integrated Spectra 

We calculate four standard model disks with a = 0.1 
accreting onto a Schwarzschild black hole with /=0.01, 
0.03, 0.1, and 0.3. Although larger values of I are as- 
trophysically interesting, we consider only models with 
I < 0.3. For I ~ 1, radial advection begins to become 
important and the assumption that all locally dissipated 
accretion power is radiated locally ceases to be accurate. 
In this regime, a slim disk model should be used in place 
of the thin disk model which is no longer self-consistent. 

In Figure 13 we plot the standard model SED's (solid 
curves) observed from infinity at i = 70°. At each value 
of I, we also plot best- fit, fully- relativistic spectra from 
Schwarzschild disks in which the local SED at each annu- 



lus is assumed to be a diluted blackbody. These diluted 
blackbody disks have the same surface brightness and 
utilize the same transfer function as the standard model 
disks. Therefore, they do not directly match the diskbb 
model which assumes a simplified form for the radial de- 
pendence of the flux and neglects relativistic transfer ef- 
fects. The best fit values of / co i = 1.4, 1.46, 1.56, and 
1.62 for the I = 0.01, 0.03, 0.1, and 0.3 models respec- 
tively. Because of limb darkening, the fully non-LTE at- 
mosphere spectra have a lower apparent luminosity than 
the isotropic diluted blackbodies at i = 70. To account 
for this, the diluted black body spectra are plotted at 
81% of their intrinsic luminosity. 

We performed the fits shown in Figure 3 by creating 
XSPEC table models from diluted blackbody spectra and 
the fully non-LTE spectrum at the same I and i. For each 
/ and i, we use the non-LTE model SED to generate an 
artificial PHA data set using a diagonal response ma- 
trix with constant effective area. This data set is then 
fit over the 0.1-10 keV range with the corresponding di- 
luted blackbody table model which has two parameters: 
/coi and the normalization. We account for the effects 
of limb darkening in the non-LTE atmospheres by fitting 
for the normalization of the diluted blackbody model and 
find that the best fit models deviate from the artificial 
spectra by less than 10-20% except out in the high en- 
ergy tail. (In an observed source, this tail would have 
low photon statistics.) If we had not fit the normaliza- 
tion the deviations would be as large as 50% at higher 
inclinations. 

Because a and M are fixed, the emitting area of the 
disk is fixed. An increase in luminosity can only result 
from an increase in the surface brightness of the disk sur- 
face. For blackbody emission, f pe ak oc T ff(r max ) cx I 1 / 4 . 
The same general trend is observed for our models, but 
the scaling does not strictly hold due to deviations from 
isotropic blackbody emission in the atmosphere calcula- 
tions. As I increases, it can be seen that standard model 
SED's become increasingly hard relative to a blackbody 
disk at the same luminosity. That is, the solid curves 
with high (low) / are better fit with diluted blackbodies 
with higher (lower) values of / co i- This relative hardening 
is largely due to the increased dominance of electron scat- 
tering in the more luminous models. The trend can be 
seen in Figure 01 where we plot the relative contributions 
of the opacity mechanisms at r* in the r = 12.6 annuli. 
As I increases, the radiative flux increases and the met- 
als become highly ionized. The corresponding decrease in 
the bound-free opacity at lower photon energies moves t* 
nearer the midplane where the temperatures are larger. 
This produces a modified blackbody spectrum which is 
harder for higher values of I. For I = 0.3, only the high- 
est ionization states of iron are sufficiently populated to 
provide significant bound-free opacity. This produces a 
larger value for r* and |/* > 1, implying that Comp- 
tonization is saturated. 

ST95 attribute the effectiveness of the diluted black- 
body approximation for their higher luminosity models 
to the presence of saturated Comptonization producing 
Wien-like spectra in the inner annuli of their disks. The 
impact of Compton scattering on our spectral calcula- 
tions can be approximately estimated from Figure [S] in 
which we plot y* in the r — 12.6 annulus of each model. 
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We find that y* > 1 only in the I = 0.3 model (solid line). 
Nevertheless, the diluted blackbody spectra still deviate 
from the non-LTE model spectra by at most 10-20% at 
typical photon energies above 0.1 keV. As discussed in 
section mi this is in part due to the relatively weak fre- 
quency dependences of T* and e v near the spectral peak. 
The agreement between the two types of spectra gener- 
ally improves with increasing I. The effects of increas- 
ingly saturated Comptonization may account for some 
of this trend. The reduction in bound-free opacity as the 
degree of ionization increases probably also plays a role 
since some of the deviations are associated with spectral 
breaks at the photoionization edges. 

In constructing these SED's we have neglected the ir- 
radiation of the disk surface by an external corona. Pro- 
vided the local ionization parameter of the irradiating 
coronal flux is sufficiently low, the reflection spectrum 
will be determined predominantly by the intrinsic ion- 
ization state of the disk. For illustration, we show the 
ionization fraction of iron at a Thompson depth of unity 
in Figure El for the I = 0.3 disk. Inside r = 100 Fe XXV 
is the dominant ion, with a small fraction of Fe XXVI 
and Fe XXVII at r > 10. The drop in ionization state 
near r ~ 6 is due to the no-torque inner boundary con- 
dition. For this disk model, we expect the Fe Ka line 
features in the reflection component to be dominated by 
the Fe XXV recombination line. Resonant trapping lim- 
its fluorescent line emission for Fe XVII through Fe XXIII 
(Ross & Fabian 1993) so any feature from more neutral 
iron must come from r > 1000. 

3.3. Dependence on Spin 

In an X-ray binary, the case of non-zero black hole spin 
seems particularly plausible. Here we consider the "max- 
imally spinning" a = 0.998 case (Thorne 1974) as an up- 
per bound on the disk spin with the expectation that 
angular momentum extraction mechanisms likely limit 
real black holes to lower spin parameters. 

Changing a modifies the one-zone radial disk structure 
equations and the location of r ms (NT73). Near the event 
horizon, the relativistic effects on the photon geodesies 
are amplified, altering the SED observed from infinity. 
At fixed I, the accretion rate drops due to the increased 
efficiency. However, the innermost annuli are still hotter 
than in the Schwarzschild case and the SED is expected 
to be harder modulo relativistic effects on the geodesies. 

We construct three maximally spinning Kerr models 
with a = 0.1 accreting at I — 0.01, 0.03, and 0.1. Diffi- 
culties obtaining convergence in some of the inner annuli 
have so far prevented us from completing an I = 0.3 
disk. In Figure [7] we plot the integrated SED's for these 
models observed at i = 70° by an observer at infinity. 
At each value of I, we also plot the best-fit, fully rela- 
tivistic spectra from maximally spinning disks in which 
the spectrum at each annulus is assumed to be a diluted 
blackbody. We fit them with the method described in 
section l3~2"l Comparison with Figure shows that both 
the full atmosphere models and the diluted blackbody 
curves are harder than their Schwarzschild counterparts 
at the same I. Increasing a modifies the one-zone ra- 
dial disk structure equations, yielding larger fluxes and 
correspondingly higher T c g in the innermost annuli. 

Relativistic broadening also accounts for the harden- 
ing of the SED. As a increases from zero to 0.998, r ms 



decreases from 6.0 to 1.23. Moving closer to the event 
horizon amplifies the relativistic effects on the photon 
geodesies, altering the SED observed from infinity. At 
i = 70° , strongly blue-shifted emission from matter mov- 
ing toward the observer pushes the high energy tail to 
larger frequencies. In Figure HO we plot the same set of 
models as viewed from a lower inclination of i — 45°. 
The relativistic broadening is clearly reduced and both 
the spectral peak and the high energy tail occur at lower 
energies. 

In the hot inner annuli, the metals arc highly ionized. 
Even for iron, a substantial fraction of the atoms are com- 
pletely stripped of electrons in the / = 0.1 model. This 
reduces the bound-free opacity and electron scattering 
plays a more dominant role than in the Schwarzschild 
case. As a result, the diluted blackbody models which 
best approximate the solid curves generally have higher 
values of / co i and the breaks near the metal edges are 
not as prominent. A calculation for these hotter annuli 
finds > 1 for most relevant photon energies in the 
I = 0.03 and 0.1 models. The agreement between the di- 
luted blackbody and atmosphere model SED's improves 
with increasing I. For i — 70°, the maximum devia- 
tions are always less than or equal to those in the cor- 
responding Schwarzschild model. Since bound-free opac- 
ity is reduced, the ratio of absorption to scattering is 
much more frequency dependent for the higher luminos- 
ity models and the improved agreement between the two 
types of spectra is likely due to Comptonization. The 
local limb darkening at characteristic photon energies 
emitted by the hot inner annuli is even better approxi- 
mated by that of a semi-infinite Thompson scattering at- 
mosphere (Chandrasekhar 1960) than the Schwarzschild 
case depicted in Figure El 

3.4. Dependence on a 

As discussed in section [21 the set of three parameters 
mo, T e ff , and Q along with the abundances and dissipa- 
tion profile, uniquely determine the structure and spec- 
trum of an annulus in our models. Of these quantities, 
only mo depends on a. If all other quantities are fixed, 
an increase (decrease) in a produces a larger (smaller) 
stress. For the accretion rate to remain constant, there 
must be a corresponding reduction (increase) in the sur- 
face density and therefore in mo- In a constant density, 
radiation pressure dominated one-zone model, a change 
in mo leads to an equivalent fractional change in the den- 
sity because the scale height is independent of a (SS73). 

In Figure we compare the vertical temperature and 
electron number density n e profiles in disks with different 
values of a. The solid and dashed curves correspond 
to a = 0.1 and 0.01 respectively. We plot results for 
annuli at r = 12.6 in disks with a = and I = 0.1. 
Both of these atmospheres are expected to be radiation 
pressure dominated so we scale the z coordinate with the 
appropriate one-zone model h (HH98 eq. 53). The x's in 
Figure|5]mark the location of r* evaluated at a frequency 
near the spectral peak (^ pC ak ~ 5 x 10 17 Hz _1 ). 

Below one scale height the models are consistent with 
the predictions of a one-zone model. The densities are 
nearly constant and n e in the a = 0.01 atmosphere is 
about a factor of ten larger than in the a = 0.1 case. The 
flux rises linearly with z providing a force that balances 
gravity. However, above z = h the flux rises less steeply 
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and begins to asymptote to crT^ s . Gas pressure gradients 
are therefore needed to balance the continued increase in 
gravity and the density decreases exponentially. The x's 
indicate that most of the spectrum is formed where the 
density is dropping rapidly and not where the density is 
relatively constant. Disk models which neglect vertical 
gradients in density will therefore overestimate the ratio 
of absorption to scattering opacity and underestimate the 
effective Compton y-parameter. 

Despite their differences for z < h, the temperature 
and density profiles seen from the surface down to t* 
are nearly identical in the two annuli. In other annuli 
at larger r and smaller T e g , the t* surfaces move further 
from the midplane so that regions of spectral formation 
remain near the surface where the density decreases ex- 
ponentially. In general, the temperature, density and 
ion level population profiles of the two models still do 
not differ significantly in these surface regions at any r. 
In FigurelTUIwe compare the integrated SED's from these 
disks viewed at i = 70°. They are nearly indistinguish- 
able. 

4. THE STRESS PRESCRIPTION 

Most of the standard model disks (a-disks) presented 
in section |21 are radiation pressure dominated in the in- 
nermost annuli. Such disks have long been known to be 
"viscously" and thermally unstable. One way to avoid 
these instabilities is to construct models (hereafter f3- 
disks) in which the accretion stress scales with gas pres- 
sure P gas rather than total pressure P in equation J2J 
(Lightman & Eardley 1974). It has been argued that 
such a scaling could be the consequence of magnetic field 
buoyancy limiting magnetic stresses (Sakimoto & Coro- 
niti 1981; Stella & Rosner 1984). 

As discussed in section EOl altering the stress prescrip- 
tion only changes m-o in a one-zone model. A process that 
limits the magnetic stress (e.g. magnetic field buoyancy) 
may also alter the vertical dissipation profile, but we dis- 
cuss the impact of the dissipation profile separately in 
section EJ For fixed a, replacing P by P gas in equation 
@ results in a lower stress. Thus, a /3-disk requires a 
larger mo to produce the same accretion rate as an a- 
disk. The larger m results in greater midplane densities 
and temperatures. 

We construct /3-disk models with a = 0.1 for 
Schwarzschild black holes with I = 0.03, 0.1, and 0.3 
and maximally spinning Kerr black holes with I = 0.03 
and 0.1. At I = 0.01 the a-disk models are gas pres- 
sure dominated down to r; n and are therefore equivalent 
to /3-disks. We calculate 7% by integrating the LTE- 
grey model until equation is satisfied for P = P gas . 
In Figure 1111 we compare the integrated SED's of the 
Schwarzschild models with the equivalent a-disk spectra 
for an observer at i = 70° . The results for the maximally 
spinning Kerr models are similar. We find that the stress 
prescription has little effect on the black hole disk SED's 
for the same reasons as in section EOl Specifically, r* lies 
near the surface of the disk where the vertical structure 
of the two types of annuli are always similar regardless 
of the conditions at the midplane. Photons emitted at 
r* in the /3-disks encounter nearly identical profiles of 
density, temperature, and level populations as those in 
an a-disk. 



5. TORQUES ON THE INNER BOUNDARY 

The no-torque boundary condition may be invalid if 
magnetic fields keep the material in the plunging re- 
gion causally connected to the disk (Gammie 1999; Kro- 
lik 1999). In that case, the surface density and bright- 
ness of the disk are modified from the standard model. 
The corresponding changes in T e g(r) and mo(r) can be 
parametrized in terms of Am the change in efficiency 
of converting energy into radiative flux (Agol & Kro- 
lik 1998). For a Schwarzschild black hole with no inner 
torque, T) = 0.057. We calculate a torqued Schwarzschild 
disk with I = 0.1, a = 0.1, and Ar/ = 0.05 so that the 
efficiency is nearly double the no-torque case. Since we 
fix the luminosity of the two disks, the accretion rate 
of the torqued disk is lower. Our method of integrat- 
ing the LTE-grey model fails to converge for some of the 
hottest annuli so we use the one-zone calculation to find 
mo throughout the disk. We computed the values of a 
required to satisfy equation[2]in the converged non-LTE 
annuli, and they differed from the assumed value by less 
than 50 percent. For the reasons discussed in section rOI 
we do not expect the the resulting spectrum to be sen- 
sitive to this small discrepancy. In the disk without a 
torque, T c g increases with decreasing r until reaching a 
maximum at r max > r ms and then falls to zero at r ms . For 
the torqued disk, r max = r ms and T c g continues to rise 
all the way down to the inner edge of the disk. Despite 
the decrease in rh, the fluxes in the innermost annuli are 
higher than in a standard disk at the same luminosity. 

In Figure ^1 we compare the SED's of the torqued 
(solid) and standard model (dotted) disks viewed at 
i = 70°. Using the same surface brightness and fold- 
ing through the same transfer function, we also calculate 
spectra assuming the local emission produces an isotropic 
diluted blackbody with a best-fit / co ] = 1.61. Owing to 
the larger fluxes, the innermost annuli of the torqued 
disk are hotter and the SED is much harder. In these 
annuli, the metals are highly ionized and provide little 
opacity. Comptonization is saturated leading to annuli 
with Wien-like spectral profiles. The diluted blackbody 
is still good to about 10% for frequencies at or below 
the peak. The torqued disk requires a slightly larger / co i 
than the corresponding standard model disk. 

One motivation for examining torqued disks is the hope 
that they might be well approximated by a diskbb model. 
Since the diskbb model is relatively inexpensive to calcu- 
late it is particularly useful for fitting data. It assumes 
a simple surface brightness profile which ignores the no- 
torque inner boundary condition and the relativistic cor- 
rections to the disk structure equations. The surface 
brightness of the torqued disk provides a better match 
to this profile so we fit our fully relativistic models with 
the diskbb model. Fits for the disk models with no inner 
torque are discussed in section 

We perform the fits by first generating XSPEC table 
models from our model SED. Next, we use these table 
models along with the XMM-Newton European Photon 
Imaging Camera pn (EPIC-pn) response matrices to cre- 
ate artificial PHA data sets. We also account for the 
intervening Galactic absorption with a neutral absorber 
model. Accounting for photoelectric absorption is im- 
portant because it affects the statistical weight of the 
bins used in the fits. The absorption reduces the relative 
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fraction of photon counts on the low energy end of the 
EPIC-pn band and makes the fit more sensitive to the 
emission at higher energies. We choose a H column of 
5 x 10 21 cm~ 2 as a representative value for a typical X- 
ray binary source. We also assume the source is at 5 kpc 
with 75 ksec exposure, providing a high signal-to-noise 
ratio dataset. We then fit these data with a diskbb model 
over the 0.3-10 keV energy range including the same ab- 
sorption used to generate the artificial spectrum. The 
reduced \ 2 of the fit is poor with ~ 20% deviations be- 
low 1 keV and ~ 10% differences near the spectral peak. 
The best-fit diskbb model is narrower than our model, 
providing a deficit of flux at the high and low energy 
ends of the band and an excess near the spectral peak. 
This is in part due to the departures from blackbody 
emission, but is due mostly to relativistic transfer effects 
and the remaining differences in the surface brightness 
profile. Considering a larger torque would increase the 
similarity in the surface brightness profiles. This should 
improve the quality of the fit but discrepancies would re- 
main. Even when we assume locally isotropic blackbody 
emission with a surface brightness identical to the diskbb 
model, the resulting SED is intrinsically broader than 
the diskbb due to the effects of relativistic transfer. 

6. THE DISSIPATION PROFILE 

In section we stated that our standard model as- 
sumes that the energy dissipation rate per unit volume e 
is locally proportional to the density p. Before exploring 
other dissipation profiles, it is worthwhile reminding the 
reader what this implies for a radiation pressure dom- 
inated disk. In that case, hydrostatic balance reduces 
to 



-1 dP r 



P 



dz 



k f F _ GM C 
c ~ ~W~B' 



(6) 



where P ra( j ~ P is the radiation pressure, C and B are 
general relativistic correction factors (Riffert & Herold 
1985), and kf is flux mean opacity. In these scattering 
dominated disks, kf k « cs which is only weakly depen- 
dent on z through the free electron fraction. Therefore 
F oc z, implying that e is constant when the gradient 
in radiation pressure balances gravity. Our assumption 
that e oc p then implies that the density is constant so 
long as gas pressure gradients are negligible. From the 
definition of column mass dm = —pdz, it is easy to see 
that our assumption reduces to 



e = p 



dF 
dm 



= P- 



mo 



(7) 



where a is the Stefan-Boltzman constant. The first 
equality holds by definition and the second follows from 
our assumption which implies dF / dm is independent of 
m. We have set F(m) — erT 4 ff (l — m/mo) in all the 
models discussed above. 

6.1. Physically Motivated Dissipation Profiles 

We would like to examine how the choice of F(m) af- 
fects the vertical structure of the accretion disk but there 
is no firm theoretical basis to guide our analysis. It is 
commonly accepted that the magnetorotational instabil- 
ity (MRI, Balbus & Hawley 1991) can provide a means 
for tapping the free energy of the differential rotation 



and producing turbulence in which energy may cascade 
irreversibly from large to small scales. However, the na- 
ture of non-radiative vertical energy transport and the 
mechanisms for dissipation in these disks are still not 
well understood. Despite these uncertainties, there is a 
rather robust prediction of magnetic field buoyancy in 
radiation dominated annuli (Sakimoto & Coroniti 1981; 
Stella & Rosner 1984). It is therefore quite reasonable 
to expect that a larger fraction of the dissipation occurs 
nearer the surface than is assumed in equation 

Numerical simulations which include radiative diffu- 
sion in a vertically stratified disk (Turner 2004) bear out 
this prediction. The radiation-MHD equations are solved 
for a patch of disk centered at 200 R g from a 10 8 M@ 
black hole. Magnetic field energy is produced fastest by 
the MRI two to three density scale heights away from the 
midplanc. The magnetic field is buoyant and rises toward 
the surface at approximately the Alfven speed. The den- 
sity is more centrally concentrated than in a standard 
model disk with the same surface density and radiative 
flux. The time averaged dissipation is not at all consis- 
tent with the e oc p assumption as much of the dissipation 
takes place in the low density surface layers. 

Similar calculations are not yet available for disks ac- 
creting onto IOMq black holes. In order to explore the 
possible consequences of magnetic buoyancy, we choose 
a form for dF / dm which approximates the Turner (2004) 
results for a 1O 8 M0 black hole. We then use this form to 
calculate F{m) with choices of too and T c g appropriate 
for an annulus around a 10M Q . Although these simula- 
tions are the best quantitative model for F(m) available, 
this choice should not be considered overly restrictive. 
There is still considerable uncertainty in the simulation 
itself. Although 29% of the dissipation is due to Silk 
damping (Agol & Krolik 1998), most is still the result 
of numerical reconnection. Furthermore, the outermost 
grid zone is quite optically thick but our models extend to 
t<1so there are several decades in to for which F{m) 
is not constrained by the simulation results. Following 
HH98, we approximate dF/dm as a broken power law 
with to being the dependent variable. The broken power 
law is parametrized by a division point to^ = /dWo and 
exponents above (Ci) and below (Co) the division. For 
the standard model, Co = Ci = 0. In the simulation, 
the dissipation and density are functions of position and 
time so we average the domain horizontally and we av- 
erage over time from 20 to 40 orbits. The broken power 
law provides a reasonable approximation with Co = —0.9, 
Ci = 0, and /d = 0.004 so we do not consider more com- 
plicated profiles. 

Next we insert this choice of dF/dm in an annulus 
with too and T e ff appropriate for r = 12.6 in a I = 0.1, 
a = 0, a = 0.1, and 10M Q standard model disk. We had 
difficulties obtaining convergence for non-LTE models in 
which the Eddington factors were updated after each so- 
lution of the radiative transfer equation. If we instead 
fix the Eddington factors after the first radiative trans- 
fer solution, the models converge. We believe the lack of 
convergence in the first case is a result of our calculation 
method rather than an instability in the model. Fur- 
thermore, we doubt the assumption in the second case 
has a strong effect on the resulting spectra. The maxi- 
mum fractional change in the Eddington factors after a 
global solution of the radiative transfer is less than 3% for 
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all frequencies and depths and is generally significantly 
smaller. 

The resulting annulus structure (dotted curve) is com- 
pared with the corresponding standard model (solid 
curve) in Figure fHfl The two are quite different. The 
radiative flux is plotted in the top panel. For z < h, the 
solid curve rises linearly because the radiative flux must 
balance the linearly increasing vertical gravity. The dot- 
ted curve falls below this line so the gradient in the gas 
pressure also provides significant support against gravity 
in the modified annulus. The temperature is shown in the 
middle panel. The surface layers above r = 1 (marked 
by the squares) are somewhat hotter in the modified an- 
nulus owing to the increased fraction of dissipation in 
that model. The electron number density is shown in 
the bottom panel. The density is more centrally concen- 
trated than in the standard model annulus as a density 
gradient is necessary to compensate for the reduced radi- 
ation pressure support. The temperatures and densities 
at r* (marked by x's) are nearly equal in the two models, 
though the shape of the temperature and density profiles 
are more dissimilar from each other than in the models 
with differing stress prescriptions. 

The specific intensities at i — 55° are plotted in Figure 
1141 The SED of the modified annulus is slightly harder 
giving an / co i which is about 10% higher and the He-like 
Fe absorption edge at 8.8 keV becomes an emission edge. 
Despite these differences, the model spectra are still qual- 
itatively quite similar. The strong temperature inversion 
in the modified annulus is above the photosphere for fre- 
quencies redward of the 8.8 keV edge and has little effect 
on the resulting SED. 

6.2. Dissipation and Coronae 

The geometry and origin of the coronae in accretion 
disks remain uncertain. One plausible explanation is 
that the coronae are due to hot electrons heated by re- 
connection of magnetic fields in the optically thin layers 
at the disk surface. This is the essence of the Svens- 
son & Zdziarski (1994) model where it is assumed that 
some fraction of released energy is transported without 
dissipation to a corona above the disk. It can be seen 
from the top panel of Figure El that verv little dissi- 
pation occurs above the Thompson photosphere in the 
models described in If). II We have considered several dif- 
ferent prescriptions for F{m) in which larger fractions of 
the dissipation occur at low m. The choice of a broken 
power law for dF/dm provides a smooth, continuous dis- 
sipation profile which is different from the Svensson & 
Zdziarski (1994) model in which dF/dm would be a step 
function. 

In Figure IT5l we show two annuli in which almost half 
of the dissipation occurs above the Thompson photo- 
sphere. The first case (solid curve) is similar to a Svens- 
son & Zdziarski (1994) model in that dF/dm is constant 
below the photosphere as in the standard model disks, 
but dF/dm is discontinuous at the photosphere. Nearly 
50% of the dissipation occurs above a Thompson depth of 
unity and the resulting annulus has a sharp rise (drop) in 
the temperature (density) just above the scattering pho- 
tosphere. In order to better approximate these gradients, 
we increase our number of depth points but memory re- 
strictions then limit us to just considering H, He, C, N, O, 
and Fe for this annulus only. The neglect of less abundant 



species seems to have little effect on the vertical structure 
or spectrum. In the second case (dotted curve) dF/dm is 
a continuous broken power law with £o=-l-l, Ci = 0> an d 
/d = 0.0001. About 40% of the dissipation occurs above 
the Thompson photosphere and another w 20% occurs 
above r* (measured near the peak). We consider the 
same mg and T c g used for the models plotted in the Fig- 
ures El an d El The profiles in the discontinuous model 
closely resemble those of a standard model disk below 
the photosphere, but above there is a sharp rise in the 
temperature and drop in the density. In the second case, 
the shape of the density and temperature profiles are 
similar to the modified annulus discussed in section f6. II 
but on a more exaggerated scale. The surface layers are 
much more extended, the surface temperature and cen- 
tral density are higher, and the densities at the surface 
are lower. 

We plot the SED's of these annuli in Figure El I n 
the top panel we plot the specific intensity at 55°. In 
both cases, Compton scattering by the hot thermal elec- 
trons above the photosphere produces a steep tail with 
emission up to ~ 100 keV. These tails are only crudely 
approximated by a power law as there is curvature at low 
energies near the peak and at high energies due to the 
thermal cutoffs. The discontinuous model (solid curve) 
produces a slightly harder spectrum. A flatter, more 
power law-like tail such as is commonly observed could 
presumably be produced by choosing a more top heavy 
F(m), but the surface temperatures become too large 
and our treatment of Compton scattering is no longer 
valid. We also include power law curves with photon in- 
dices of r = 3.4 and 3.5 to approximate the two spectra 
over the ^5 — 50 keV band. In the bottom panel we 
plot the ratio of the model SED's to these power laws. 
For the discontinuous model there is curvature due to 
the reduction in the scattering albedo from bound-free 
opacity above ~ 10 keV and the spectrum resembles the 
putative Compton reflection hump. The high energy side 
of the hump is due to both Compton down scattering by 
the cooler electrons near the photosphere and the ther- 
mal cutoff of the hot electrons. There is a strong emis- 
sion feature at the He-like Fe edge. This feature is likely 
enhanced because our treatment only allows recombina- 
tions to the ground state. Some fraction of this power 
would be shifted to lower energies if we had considered 
recombinations to excited levels. In the second case, the 
spectrum simply falls off at high energies and there is no 
structure in the residuals which could be characterized 
as a reflection hump. A large fraction of dissipation oc- 
curs in the spectral forming regions just below the photo- 
sphere so the temperature there is high and the density 
is relatively low. Almost all iron is completely ionized 
so that the bound-free opacity which produces the low 
energy side of the hump in the previous model is insignif- 
icant. 

7. LUMINOSITY - TEMPERATURE RELATION 

Detailed fits to individual observations are beyond the 
scope of this paper. Because we neglect irradiation of 
the disk surface, we focus on comparison with the results 
of GD04. Their sample is particularly relevant because 
they select observations in which the coronal emission is 
less than 15% of the bolometric luminosity of the source. 

GD04 construct plots of the temperature-luminosity 
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relation for each source in their sample. The luminosity 
of the disk model Ldisk is plotted as a fraction of £Edd 
versus the maximum color temperature T max . A major 
difficulty with inferring the physical properties of accre- 
tion disks is that the distance, inclination and black hole 
mass of the source are often not all known with great 
accuracy. Therefore, the vertical and horizontal position 
of the locus of points on these curves is uncertain. How- 
ever, the shape of the locus of points for an individual 
source should be robust to such uncertainties. 

If the disk spectrum is well approximated by a mul- 
titemperature diluted blackbody with fixed f co \, Ldisk 
should be proportional to T^ ax . Specifically Gierlihski 
et al. (1999) assume that the gravitational field may 
be approximated with a pseudo-Newtonian potential 
(Paczyhski & Witta 1980) and find the relation 




where T max = T e s imax / co i ■ GD04 find that most of the 
sources roughly follow this scaling, though several show 
some degree of relative hardening (increasing / co i) as 
Ldisk increases. 

One obstacle to comparing with their results is that 
neither T max nor / co i is a well defined theoretical quan- 
tity in our models. In principle, we could use the / co i val- 
ues from the best- fit models shown in Figure |3 but that 
fitting procedure is quite different from the method im- 
plemented by GD04. In constructing their temperature- 
luminosity diagram, they fit the diskbb model to their 
spectra to account for the presumed disk component. 
They then use the best fit inner temperature and to- 
tal model flux to calculate idisk/^Edd and T max . In do- 
ing so they apply relativistic correction factors (Zhang, 
Cui & Chen 1997) for a Schwarzschild black hole at the 
appropriate inclination. They also apply an additional 
correction factor to account for the mismatch between 
the surface brightness profile of the diskbb model and the 
pseudo-Newtonian potential used to derive equation [SJ 

To facilitate comparison, we try to reproduce this pro- 
cedure as closely as possible. We produce an artificial 
EPIC-pn PHA data set as described in section [5] We 
also generate a second set of artificial spectra by the same 
method but employing the RXTE Proportional Counter 
Array (PCA) response matrices. Galactic absorption has 
less impact on the PCA photon statistics so we do not in- 
clude it in those artificial spectra and fits. The EPIC-pn 
spectra are useful because of the high signal-to-noise and 
the low energy coverage while the PCA spectra are neces- 
sary for comparing with GD04. Next, we fit diskbb mod- 
els over the 0.3-10 keV and 3-20 keV energy ranges for the 
EPIC-pn and PCA data sets respectively. We also gen- 
erate artificial spectra including a compTT (Titarchuk 
1994) component which is 10% of the bolometric flux 
to account for possible Comptonized coronal emission. 
These spectra are then fit with a combined diskbb + 
compTT model. We find that the best fit parameters of 
the diskbb models are only weakly sensitive to the pres- 
ence or absence of the compTT component as long as we 
fit for this component with another compTT model and 
not a power law. Finally, we apply the same prescription 
outlined in GD04 to calculate idisk/^Edd and T max from 
the model temperature and flux. 



We show the results of this procedure for the four 
Schwarzschild models with I = 0.01, 0.03, 0.1, and 0.3 
viewed at i — 45° and 70° in Figure El The solid curves 
represent the luminosity-temperature relation of equa- 
tionElfor f mX = 1.4, 1.6, 1.8 and 2.0. The triangles and 
squares mark the XMM-Newton EPIC-pn and the RXTE 
PCA measurements respectively. Only the highest three 
luminosities are plotted for the PCA measurements be- 
cause the model with / = 0.01 emits only a small fraction 
of its luminosity above 3 keV and is not well constrained 
by a PCA observation. The diskbb model provides a poor 
fit to the artificial XMM data so reliable uncertainties 
cannot be estimated for best fit parameters. The un- 
certainties (90% confidence for one parameter) on PCA 
depend on the signal-to-noise of the simulated data set, 
but they are of order the symbol size or smaller for a 1.5 
ksec exposure of a source at 5 kpc. We compare the best- 
fit diskbb spectra to the i = 70° non-LTE model SED's 
in Figure ITSl 

First, we focus on the EPIC-pn data sets. The filled 
and open triangles mark the i = 45° and 70° measure- 
ments respectively. Both sets of measurements show a 
similar degree of spectral hardening with increasing lu- 
minosity but with approximately the same fractional off- 
set in T max and idisk/^Edd for each I. The difference 
in apparent luminosity between the two sets of measure- 
ments and the known I values are predominantly due to 
limb darkening in the atmospheres (compare with Fig- 
ure [21, but there is also a small mismatch between our 
relativistic transfer function and the g correction factors 
we interpolate from Zhang et al. (1997) Table 1. Com- 
paring the symbols with the solid lines shows that the 
derived values of / co i are inclination dependent. When 
folded through the spectral response, the best-fit diskbb 
models disagree with artificial spectra by about 5-10% 
at lower energy and less than 5% near the peak. These 
deviations can be significant for the high signal-to-noise 
obtained in a long exposure of the bright high/soft state 
X-ray binaries. The model provides poor fits (reduced 
X 2 S> 1) to the data for the simulated 75 ksec exposures 
of sources at 5 kpc. 

The PCA measurements at I = 0.1 and 0.3 agree rea- 
sonably well with the EPIC-pn data sets. The T max 
values are larger because the PCA band only extends 
down to 3 keV in our fits and does not include the spec- 
tral peak. The Doppler-broadened high energy tail is 
consistent with a slightly higher inner temperature than 
the emission near the spectral peak which is given more 
weight in the fits to the EPIC-pn artificial spectra. The 
I = 0.03 case is slightly more complicated because there is 
a spectral break above 3 keV due the presence of bound- 
free opacity. This produces a slightly steeper fall-off in 
the high energy tail. This can be more easily fit by a 
lower temperature and correspondingly higher normal- 
ization leading to a larger idisk/^Edd and lower T max . 
These shifts overestimate the emission near the peak 
which is outside of the PCA band but must be accounted 
for in the EPIC-pn fits. 

8. DISCUSSION 

8.1. Comparison with Previous Work 

Our models are qualitatively different from those of 
ST95 in several important ways. We examine the specific 
intensities and account for relativistic effects on photon 
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geodesies while they compare the specific flux in the local 
frame of the disk. Calculating specific fluxes in the local 
frame simplifies the comparison with diluted blackbodies 
by avoiding inclination-dependent effects. However, this 
has the drawback of neglecting the non-trivial coupling 
between frequency and inclination caused by the limb 
darkening and relativistic transfer which is presumably 
present in real accretion disks. Another important differ- 
ence between our models is our inclusion of metals. We 
account for the bound-free opacity of all important ions 
in addition to free-free opacity which was the only ab- 
sorption opacity included in their calculations. The addi- 
tion of metal opacity reduces the range of I in which sat- 
urated Comptonization is important and produces modi- 
fied blackbodies which are slightly softer (more consistent 
with lower values of / co i) than their model SED's. 

ST95 conclude that multi-temperature, diluted black- 
body spectra adequately approximate their models at 
higher luminosities but are inadequate at low luminos- 
ity. When they compare diluted blackbodies with their 
a = 0.1, Schwarzschild black hole integrated disk spec- 
tra at energies above 0.1 keV, the deviations are less than 
20%, 15%, and 30% for the I = 0.0057, 0.057, and 0.57 
models respectively. For the two more luminous models, 
the largest disagreement is at the low energy end of the 
band and the agreement is better at higher energies. For 
I = 0.0057, the integrated spectra show larger deviations 
at higher energies and the individual model annuli are 
not well represented by diluted black bodies because the 
y-parameter is small. They conclude that color-corrected 
blackbodies are good approximations only at the higher 
luminosities. 

We focus on the integrated emission from the whole 
disk. For all the I = 0.01 to 0.3 models, the com- 
parisons between diluted blackbody and non-LTE atmo- 
sphere model SED's depend on the viewing angle. In sec- 
tion |^21 we presented fits to our non-LTE models with 
fully relativistic diluted blackbody models with constant 
color correction f co \. Provided one fits for both f co \ and 
the normalization (in order to account for limb darken- 
ing), these fits are surprisingly good. The disagreement 
between the models is less than 20% at most photon 
energies regardless of inclination. As seen in Figure |3| 
the deviations in the high energy tails can eventually be- 
come much larger. However, by that point the statistical 
weight of these energy bins is typically low because the 
integrated photon flux falls off rapidly in the tail. The 
fitted values of f co \ at i = 70° are summarized in the sec- 
ond column of Table These best-fit values are weakly 
dependent on inclination, but they are always lower than 
the 1.7-1.9 values reported by ST95. 

A direct comparison of our models with those of Mer- 
loni et al. (2000) is complicated by their assumption 
that a fraction / of the dissipation occurs in the corona. 
Inspection of their Figure 2 suggests that the value of 
/coi inferred from their spectral models is mostly deter- 
mined by the corrected disk luminosity parameterized by 
(1 — f)m. Specifically, models with different / but the 
same (1 — f)rh yield similar values for / co i. Their m 
accounts for the disk efficiency so this disk luminosity 
parameter is the rough equivalent of our I. Merloni et al. 
(2000) also construct artificial spectra from their model 
SED's using an RXTE response matrix. They fit these 
spectra with diskbb models and derive values of f co \ by 



a method which is similar, but differs in detail from the 
procedure implemented by us in section Hand in GD04. 
They find a trend toward increasing / co i with decreas- 
ing (1 — f)rh. If we take the diskbb fits of section and 
use their method for determining / co i, we find the values 
listed in the fourth column of Tabled (For completeness, 
we also show values appropriate for the XMM-Newton 
EPIC-pn camera in the third column.) We find a weak 
increase in / co i with increasing I which is inconsistent 
with their results. We suspect that this discrepancy is 
largely due to their neglect of bound-free metal opacity 
which causes them to overestimate the ratio of scatter- 
ing to absorption opacity. Therefore t v is lower and r* 
is greater, providing a higher T*. The higher T* and 
lower e v would both give rise to harder, more modified 
spectra, consistent with larger values of / co j. These ef- 
fects are most significant in the low accretion rate models 
where Comptonization is negligible at most radii. Their 
assumption of constant density likely overestimates the 
density in the spectral forming regions in both gas pres- 
sured and radiation pressure dominated annuli, increas- 
ing the ratio of absorption to scattering opacity at r*. 
This reduces the y-parameter but may somewhat miti- 
gate the increased spectral hardening which comes from 
neglecting metal opacity. 

8.2. Stress Prescription and Dissipation Profile 

For the range of I, a, and a explored here, altering 
mo by varying a or changing the stress prescription has 
little effect on the SED. Even changes in the vertical 
dissipation profile do not matter very much unless a sub- 
stantial fraction of the heating occurs above the effective 
photosphere. It therefore appears that the spectra are 
insensitive to the details of the very uncertain vertical 
structure of disks, which might be encouraging from a 
theoretical point of view, but disappointing if one wants 
to use spectral observations to probe disk physics. 

However, it is important to bear in mind how this ro- 
bustness arises. The innermost parts of luminous black 
hole accretion disks are both radiation pressure sup- 
ported and electron scattering dominated. Provided heat 
is vertically transported by radiative diffusion, then the 
scale height h, which determines the surface gravity g, 
is completely independent of the stress prescription. If 
there is little dissipation in the surface layers, the radia- 
tive flux is constant, but the gravity still increases with 
height. As a result, gradients in gas pressure are always 
important in maintaining hydrostatic equilibrium near 
the surface, and this always produces a very steeply de- 
clining density profile there. This is in contrast to the 
temperature profile, which is relatively constant in the 
surface layers (cf. Fig. |5J), and close to the effective 
temperature. The steep decline in density occurs over a 
gas pressure scale height H g ~ P gas /(pg) (Hubeny 1990), 
which is independent of the details of the overall vertical 
structure as it only depends on the surface temperature 
and the local gravity. The effective photosphere at char- 
acteristic frequencies is determined by 

r cfr = J [3K th ( K ° s + K th )}^ 2 pdz = 1. (9) 

Because K th oc p and n cs oc p°, this integral is approx- 
imately proportional to p* 3 ^ 2 H g . Hence the density p* 
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at the effective photosphere is approximately indepen- 
dent of the details of the vertical structure. The col- 
umn mass down to the effective photosphere m* ~ P*H g 
and the temperature at the effective photosphere T* ~ 
T e ff(«; es m*) 1 / 4 will also be approximately independent of 
the vertical structure. The resulting spectrum will there- 
fore be quite robust. This all fails, however, if heat is 
transported by other means, such as convection, rather 
than radiative diffusion. It may also fail if substantial 
density inhomogeneities are present within the disk tur- 
bulence, as we discuss below. 

The robustness of the SED's to changes in the stress 
is also a feature of accretion disk models around super- 
massive black holes, and for the same reason. Figure 17 
of Hubeny et al. (2000) show that accretion disk models 
around a 5 x 10 s M Q Kerr hole are very insensitive to 
the value of a for Eddington ratios up to at least I ~ 0.1. 
The only discrepancies occur above the He II edge. Sim- 
ilar differences occur in stellar mass black hole models 
near the iron K edge, but this lies in the high energy tail 
of the spectrum. For high black hole spins and higher 
Eddington ratios, the effective optical depth of the in- 
nermost annuli is reduced, and the effective photosphere 
at all frequencies moves deeper into the disk atmosphere. 
In such cases the spectrum will start to depend on the 
stress and stress prescription, with high values of a or 
an a-disk prescription producing harder spectra than low 
values of a or a /3-disk prescription. The dependence on 
a is illustrated for I — 0.3 disks accreting onto super- 
massive Kerr holes in Figure 16 of HBKA. We therefore 
expect that high Eddington ratio disk models around 
stellar mass Kerr holes to be less robust to changes in 
the stress prescription as well, but we have not yet suc- 
ceeded in constructing such models. 

The weak dependence of the model SED's on the choice 
of stress prescription casts some doubt on GD04's spec- 
ulation that their luminosity-temperature relations are 
inconsistent with /3-disk models. Their reasoning relies 
on the assumption that Comptonization is not important 
in /3-disks due to the larger midplane densities inferred 
from a one zone model (Nannurelli & Stella 1989). How- 
ever, the effective photospheres in our models always oc- 
cur in the surface layers where the densities are much less 
than the midplane density. Therefore, Comptonization 
is equally important (or unimportant depending on I) for 
both a-disks and /3-disks over the range of parameters we 
consider. 

In an attempt to inject first principles physics that 
is not present in the standard model, we also consid- 
ered modifications to the vertical dissipation profile. The 
standard model dissipation profile produces a disk which 
is hydrodynamically unstable to convection in the ra- 
diation pressure dominated regions (Bisnovatyi-Kogan 
& Blinnikov 1977, HBKA), though it is far from clear 
whether significant convective transport of heat would 
occur in the presence of MRI turbulence. Magnetic field 
buoyancy suggests that the dissipation may be more con- 
centrated near the surface than in the standard model. 
An annulus with a modified dissipation profile motivated 
by simulations of accretion disk turbulence (Turner 2004) 
produces a vertical structure that is increasingly gas pres- 
sure supported in its interior, and hydrodynamically sta- 
ble to convection. 



Adopting this modified dissipation profile only pro- 
duces slight changes to the resulting spectrum, however. 
The reason is that the considerations of the first para- 
graph of this section still apply. The surface density pro- 
file is still very steep, and the physical properties of the 
plasma at the effective photosphere are largely indepen- 
dent of the dissipation profile, provided most of that dis- 
sipation occurs deeper than the effective photosphere. 
The impact is not completely negligible as the modi- 
fied annulus SED is slightly harder than in the standard 
model, though the inferred / co i only differs by ~ 10%. 
However, if we consider more top heavy dissipation pro- 
files in which a significant fraction of the dissipation oc- 
curs near or above the photosphere we can produce atmo- 
spheres with very different SED's even though the surface 
radiative flux, gravity, and surface density are the same. 
The "corona" -producing annuli discussed in section IIOl 
are extreme examples. These calculations may allow for 
a self-consistent treatment of the hard and soft spectral 
components in models in which the coronal emission is 
provided by flaring regions above the disk though there 
are several caveats to our treatment which remain to be 
addressed. These include the assumptions that the tem- 
peratures of the ions and electrons are equivalent, the 
electrons are entirely thermal, and the corona is time- 
steady. 

Most of the uncertainties in the vertical distribution 
and nature of the dissipation are relatively easy to ad- 
dress with the methods presented here. However, other 
potentially important aspects of the dynamic nature of 
these disks remain unaccounted for in our models. Sim- 
ulations (e.g. Turner 2004) suggest the presence of sub- 
stantial inhomogeneities in density due to the compress- 
ible nature of the turbulence (Turner, Stone & Sano 
2003) and radiatively driven instabilities (Gammie 1998). 
Modeling of the emergent spectrum from such structures 
requires three dimensional radiative transfer. We have 
reported on Monte Carlo radiative transfer calculations 
on inhomogeneous structures created by simulations else- 
where (Davis et al. 2003). The temperature regime was 
more appropriate for a quasar disk rather than an X- 
ray binary, and thermal Comptonization was negligible. 
The main effect of the inhomogeneities was to reduce the 
modified blackbody effects and produce a softer, more 
thermal spectrum. The emergent flux through an inho- 
mogeneous medium is also larger than the spatially aver- 
aged flux required to support the medium against grav- 
ity (Begelman 2001). This may reduce the vertical scale 
height of the disk compared to the standard model, and 
make the emergent spectrum more sensitive to the stress 
prescription and dissipation profile. Because the turbu- 
lence is compressible and supersonic, bulk Comptoniza- 
tion can also significantly affect the emergent spectrum 
(Socrates, Davis & Blaes 2004). In addition to the effects 
of Silk damping, bulk Comptonization may also provide 
a radiative dissipation mechanism for the turbulence it- 
self. It is also possible that the time-averaged spectrum 
emerging from a time-dependent, turbulent annulus dif- 
fers from the spectrum emerging from a time-average of 
the vertical structure. 

8.3. Observational Implications 

Our luminosity-temperature relation agrees reasonably 
well with the results of GD04. For example, between 
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I = 0.3 and 0.1 we find weak evolution toward lower / co i 
consistent with that observed in XTE J1550-564. The 
larger decrease in / co j in our models from I = 0.1 to 0.03 
is not clearly observed in their data, but it is not in- 
consistent due to large uncertainties in the fits at lower 
luminosities. The magnitudes of / co i derived using the 
procedure of GD04 depend on the inclination due to limb 
darkening. In this case, they are in approximate agree- 
ment with the observed values for i = 70, which is near 
the measured inclination of 73.5±^ for XTE J1550-564 
(Orosz et al. 2002). The values of / co i are probably also 
sensitive to our assumptions about inner torques, dis- 
sipation profiles, and metal abundances. Despite these 
concerns, we expect the generally weak trend toward de- 
creasing /coi with decreasing I to be robust because it 
is largely due to the increasing ratio of absorption to 
scattering opacity as the temperature of the disk drops. 
This reduces e„ and moves r* closer to the photosphere, 
lowering T*. The spectra soften and the associated / co i 
values decrease. 

Though we employ the concept of a color correction 
throughout this paper for the purposes of comparison, 
we caution that isotropic, diluted blackbodies are not 
perfect approximations to our models. Given the high 
signal-to-noise spectra afforded by an observation with 
XMM-Newton, the < 20% discrepancies between the best 
fit, fully relativistic diluted blackbodies and the non-LTE 
atmosphere spectra are potentially significant. Any as- 
sumption of isotropy is clearly poor as limb darkening 
can alter the inferred luminosity of a source by up to 50% 
at high inclinations. We further caution that the diskbb 
model is an even poorer match to our model SED's. For 
high signal-to-noise observations we do not expect it to 
provide adequate fits unless the disk is truncated at large 
radius. Even if the presence of an inner torque makes 
the surface brightness profile more consistent with the 
diskbb profile, the model completely neglects relativistic 
transfer effects. The best-fit diskbb spectra are always 
significantly narrower than our models due to differences 
in the surface brightness profiles and relativistic broad- 
ening. To accurately model all the emission from a disk 
component, one needs to include the relativistic effects 
on the disk structure and transfer function in the spectral 
model before fitting it to the observed SED. The data in 
Table ^ indicate that different but seemingly consistent 
methods of calculating / cc j with diskbb model fits give 
differing results. Furthermore, the results depend on the 
detector response and band. 

These concerns are particularly relevant for higher 
signal-to-noise observations which have inferred relativis- 
tically broadened Fe Ka lines but fit the disk compo- 
nent with an unbroadened diskbb model (Miller et al. 
2004a, 2004b). In some cases the iron line profiles are 
so broad and asymmetric that they require a near max- 
imally spinning black hole with a very steep iron line 
radial emissivity profile. Contrary to our expectations 
of 10-15% discrepancies, the diskbb model seems to ac- 
count surprisingly well for the presumed disk component 
in these observations. We suspect this may be due in 
part to their use of a power law to model the coronal 
emission. If the disk provides the seed photons for the 
corona, a power law will likely overestimate the flux of 
this component at the low energy end of the band where 



we predict the discrepancies are largest. This extra emis- 
sion could account for the difference between our model 
SED's and the diskbb models. We test this hypothesis 
by constructing an artificial spectrum with our I = 0.1, 
a = 0.998, and i = 45° model SED, assuming it is located 
at 5 kpc and accounting for an intervening H column of 
5 x 10 21 cm~ 2 . We also include a compTT component 
whose flux is about 10% of the model flux. First, we fit 
these data with a model that includes a diskbb compo- 
nent, a compTT component with the seed temperature 
fixed at the input value, and neutral absorption at the in- 
put H column. We get a poor fit with residuals of ~ 15%. 
When we replace the compTT model with a power law, 
the residuals of the best fit are reduced to less than 5% 
and the quality of fit improves dramatically. It there- 
fore seems plausible that the combinations of a diskbb 
and power law may effectively reproduce the SED of a 
real disk even though the diskbb model does not properly 
account for the accretion disk emission. 

9. CONCLUSIONS 

We have calculated several representative non-LTE, 
fully relativistic models for geometrically thin disks ac- 
creting onto a 10M© black hole. We first explored the be- 
havior of standard model (SS73, NT73) accretion disks. 
We found that the inclusion of abundant metals has a 
significant impact on the model SED's because of bound- 
free absorption opacity at high photon energies. The in- 
clusion of this additional opacity generally decreases the 
importance of Comptonization, but still produces softer, 
more Wien-like spectra than would be expected from 
free-free opacity alone. The associated diluted blackbody 
color corrections are generally lower than found by pre- 
vious authors (ST95, Merloni et al. 2000). 

The frequency dependence of the disk integrated model 
SED's is qualitatively similar to that of spectra produced 
by assuming isotropic diluted blackbodies with constant 
color correction for the local emission, provided these lo- 
cal spectra are folded through a relativistic transfer func- 
tion. Quantitative differences still exist which should be 
discernible for some observations with modern X-ray ob- 
servatories. The effects of limb darkening are significant 
and the assumption of isotropy of the local emission must 
be compensated by fitting the normalization as well as 
the color correction. The best-fit / co i's are weakly depen- 
dent on the black hole spin and accretion rate, generally 
increasing as either parameter increases. 

We have found that the standard model SED's are 
strongly affected by changes in accretion rate and black 
hole spin, but they are only weakly dependent on the 
value of a for the parameter space explored here. Modi- 
fying the stress prescription by replacing P with P gas in 
equation |21 also has little effect on the model SED's, con- 
trary to the supposition of GD04. The reason for this is 
that the overall disk scale height in the inner, radiation 
pressure supported regions is largely independent of the 
stress prescription, and so are the resulting physical con- 
ditions at the effective photosphere. Higher Eddington 
ratio models, particularly those around maximally spin- 
ning holes, may be more sensitive to the treatment of the 
stress. Steep density gradients generally exist in the sur- 
face layers, and modeling the disk annuli with constant 
density slabs is also a poor approximation. 

We found that our model SED's are generally not well 
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matched by a simple diskbb model over the whole 0.3-10 
keV band. Fits using the XMM-Newton EPIC-pn re- 
sponse have discrepancies of greater than 10%. Such 
large discrepancies are generally not seen in diskbb fits to 
high signal-to-noise observations of X-ray binaries with 
relativistically broadened Fe Ka lines. We speculate that 
the lack of residuals in these fits may be due to overesti- 
mation of the coronal emission at low energies. A resolu- 
tion of this issue will require detailed fits of our models 
to real observations. 

We used the method of GD04 to construct a 
luminosity-temperature relation for our Schwarzschild 
disk models. The luminosity scales as the fourth power of 
the temperature, but we also found a trend toward larger 
,/coi with increasing / which is generally consistent with 
their analysis for several sources. The agreement with 
their measured values of / co i is somewhat more uncertain 
and our / co i measurements depend on the inclination due 
to limb darkening effects. The values of / co i which we in- 
fer using GD04's method are higher than those we derive 
from fits with the relativistic diluted blackbody models 
(see Table 1). 

We have also looked at modifications of the stan- 
dard model based on recent theoretical progress in un- 
derstanding magnetohydrodynamic stresses in the flow. 
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TABLE 1 
Color Corrections 



/fr 



J MFR. 



fR 

J MFR 





0.3 
0.1 
0.03 
0.01 



1.64 
1.56 
1.46 
1.40 



1.62 
1.60 
1.50 
1.29 



1.67 
1.63 
1.36 



1.88 
1.85 
1.74 
1.50 



1.97 
1.90 
1.58 



Note. — Calculations of / co i based on fits to our model spectra for a = 0.1 disks around Schwarzschild black holes, viewed at an 
inclination of 70°. Each column uses a different method to determine / co i (see text). The second column summarizes / co i from our fits to 
the non-LTE model SED's with our fully relativistic diluted blackbody model spectra. A superscript X or R indicates the / co i values were 
calculated using the diskbb fits to the XMM-Newton EPIC-pn or RXTE PCA data sets respectively. A subscript MFR or GD04 means 
that we use the methods described in Merloni et al. (2000) or GD04 (respectively) to calculate / co i. 
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Fig. 1. — Specific intensity and other frequency dependent parameters evaluated in the local frame at r = 12.6 for a disk with a = 0.1, 
I = 0.1, and a = 0. (a) We compare the specific intensity viewed from 55° (solid curve) with a diluted blackbody at the same T e fj and 
/col = 1-56 (dotted curve). The units of the ordinate are ergs s _1 cm -2 ster - 1 . (b) We show the fraction of the total opacity provided by 
electron scattering (solid curve), free-free absorption (dotted curve), and bound-free absorption (dashed curve) evaluated at m*. (c) The 
depth of formation (eq. |1J is plotted, (d) The frequency dependent {/-parameter (eq. |SJ is plotted. 
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Fig. 2. — Normalized specific intensity of the emission at at r = 12.6 in an I = 0.1, c? = 0.1, and a = disk model. The curves are 
normalized by a factor of F v /n so that isotropic emission would have an ordinate of unity. The curves are plotted for photon energies 
of 0.21, (dashed) 2.2. (dotted), and 8.4 keV (solid). For comparison we plot the limb darkening law (x's) for a semi-infinite Thompson 
scattering atmosphere assuming a Rayleigh phase function (Chandrasckhar 1960). 
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Fig. 3. — Integrated disk SED's of four values of I (solid) viewed at i = 70° by an observer at infinity. At each I, we also plot the 
best-fit, fully-relativistic disk model spectra in which the local flux is assumed to be an isotropic diluted blackbody. Each model has 
a = 0.1 and a = and from lower left to upper right the curves correspond to I = 0.01, 0.03, 0.1, and 0.3. with best fit / co i = 1.4 (triple 
dot-dashed), 1.46 (dot dashed), 1.56 (dashed), and 1.64 (dotted) respectively. The diluted blackbody spectra are plotted at 81% of their 
intrinsic luminosity to account for limb darkening in the non-LTE atmosphere models. The lower panel shows the ratio of the non-LTE 
model spectra to the diluted blackbody spectra. 
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Fig. 4. — Fraction of the total opacity provided by electron scattering (solid), free-free absorption (dotted), and bound-free absorption 
(dashed) evaluated at m* and r = 12.6 in an a = 0.1 and a = disk. From top to bottom, the panels correspond to i=0.3, 0.1, 0.03, and 
0.01. 
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Fig. 5. — Frequency dependent {/-parameter, y*, for models with four different values of / evaluated in the local frame of the disk at 
r=12.6.. The curves correspond to Z=0.3 (solid), 0.1 (dotted), 0.03 (dashed), and 0.01 (dot-dashed). 
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Fig. 6. — Ion fraction of iron at the Thompson photosphere as a function of radius for an I = 0.3 and a = 0.1 disk accreting onto a 
Schwarzschild black hole. The curves correspond to Fe XXVII (solid), Fe XXVI (dotted), Fe XXV (short dashed), Fe XXIV (dot-dashed), 
Fe XVII - Fe XXIII (triple dot-dashed), and Fe I - Fe XVI (long dashed). 
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Fig. 7. — Integrated disk SED's for three values of I (solid) observed from an i = 70° by an observer at infinity. At each I, we also plot 
the best-fit, fully relativistic disk models in which the local flux is assumed to be an isotropic diluted blackbody (dashed). Each model has 
a = 0.1 and a = 0.998 and from lower left to upper right the curves correspond to I = 0.01, 0.03, and 0.1 with best fit / co i = 1.47, 1.53, 
and 1.6 respectively. The diluted blackbody spectra are plotted at 95% of their intrinsic luminosity to account for limb darkening in the 
non-LTE atmosphere models. 
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Fig. 8. — Integrated disk SED's for three values of I (solid) observed from an i = 45° by an observer at infinity. At each I, we also plot 
the best-fit, fully relativistic disk models in which the local flux is assumed to be an isotropic diluted blackbody (dashed). Each model has 
a = 0.1 and a = 0.998 and from lower left to upper right the curves correspond to I = 0.01, 0.03, and 0.1 with best fit / co i = 1.4, 1.41, and 
1.47 respectively. The diluted blackbody spectra are plotted at 114% (i = 0.01) and 112% (Z=0.03 and 0.1) of their intrinsic luminosity to 
account for limb darkening in the non-LTE atmosphere models. 
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Fig. 9. — Temperature (top) and electron number density (bottom) as a function of height above the midplane. The atmospheres arc 
located at r = 12.6 in a disk accreting on to a Schwarzschild black hole with I = 0.1 and a = 0.1 (solid) or 0.01 (dotted). The x's mark the 
position of r* evaluated for photon frequencies near the peak (^ pea k ~ 5 X 10 17 Hz -1 ) of the SED. The scale height, h, is evaluated using 
a radiation pressure dominated one-zone calculation (HH98 eq. 53). 



Disk Models of Black Hole X-ray Binaries 



23 



30 
o 



37 



CO 



bu 36 
o 




0.1 



1.0 

Energy (keV) 



10.0 



Fig. 10. — Integrated SED's observed at i = 70° by an observer at infinity for disk with I = 0.1 and a = 0. The curves with a = 0.1 
(solid) and 0.01 (dashed) are nearly indistinguishable. 




Fig. 11. — Integrated disk SED's for three values of I observed at i = 70° by an observer at infinity. Each model has a = 0.1 and a = 
and from lower left to upper right the curves correspond to I = 0.03, 0.1, and 0.3. The Q-disks (solid) are equivalent to the curves plotted 
in Figure |3] The /3-disks (dashed) show little deviation from the o-disks. 
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Fig. 12. — SED (solid) from a torqued disk with I = 0.1, a = 0, a = 0.1, and A?? = 0.05 viewed at i = 70° by an observer at infinity. We 
also plot the best-fit, fully relativistic disk models in which the local flux is assumed to be an isotropic diluted blackbody with / co i = 1.61 
(dashed). The diluted blackbody spectra are folded through the same transfer function and viewed at the same inclination as the solid 
curve. The SED of the torqued disk is considerably harder than the SED of of disk without a torque (dotted). 
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Fig. 13. — Radiative flux (top), temperature (middle) and electron number density (bottom) as a function of height above the midplanc. 
The atmospheres have Q, mo, and T c g appropriate for r = 12.6 in a disk accreting on to a Schwarzschild black hole with a = 0.1. We 
compare the standard model (Co=Cl=0; solid curve) to a model with a modified dissipation profile (£o=-0.9, Cl=0,/d = 0.004; dotted 
curve). The squares mark the position of the Thompson photosphere and the x's mark t* evaluated for photon frequencies near the peak 
(fpeak ~ 5 X 10 17 Hz _1 ) of the SED. The scale height, h, is evaluated using a radiation pressure dominated one-zone calculation (HH98 
eq. 53). 
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Fig. 14. — Specific intensity viewed at i = 55° in the local frame for the annuli shown in Figure UHl The models parameters are fo=fi=0 
(solid curve), fo=-0.9, £i=0, and f d = 0.004 (dotted curve). 
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Fig. 15. — Radiative flux (top), temperature (middle) and electron number density (bottom) as a function of height above the midplanc. 
The atmospheres have Q, mo, and T c ff appropriate for r = 12.6 in a disk accreting on to a Schwarzschild black hole with a = 0.1. 
The first model (solid curve) is similar to a standard model annulus (£o = £i = 0) except that dF/dm is not continuous and 50% of 
the flux is dissipated above the Thompson photosphere. The second model (dashed curve) has a continuous dF/dm with £o=-l.l, Ci=0, 
and /<j = 0.0001. About 40% of the dissipation occurs above Thompson depth unity. The squares mark the position of the Thompson 
photosphere and the x's mark t* evaluated for photon frequencies near the peak (^ pea k ~ 5 X 10 17 Hz -1 ) of the SED. The scale height, h, 
is evaluated using a radiation pressure dominated one- zone calculation (HH98 eq. 53). 
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Fig. 16. — Specific intensity viewed at i = 55° in the local frame for the annuli shown in Figure Il5l The first model (solid curve) 
is similar to a standard model annulus (£o = Ci = 0) except that dF/dm is not continuous and 50% of the flux is dissipated above the 
Thompson photosphere. The second model (dashed curve) has a continuous dF/dm with Qq=-1.1 } £i=0, and /<j = 0.0001 
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Fig. 17. — Luminosity-temperature relation measured for four Schwarzschild disks viewed at i = 45° (filled symbols) and 70° (open 
symbols). The solid curves represent the lines of constant f co \ = 1.4, 1.6, 1.8, and 2.0 (from top to bottom) from equation 181 The triangles 
and squares mark the XMM EPIC-PN and RXTE PCA measurements respectively. 
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Fig. 18. — Integrated non-LTE disk SED's from Figure l3l (solid curve). At each I, we also plot the best-fit diskbb models for the 
XMM-Newton EPIC-pn (dashed) and RXTE PCA (dotted) artificial data sets. The vertical solid and dashed lines mark the endpoints of 
EPIC-pn and PCA energy ranges respectively. 



